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The coexistence state of antiferromagnetism (AF) and superconductivity (SC) has been 
observed in five-layered cuprates. However, this coexistence state disappears, and the AF 
phase and SC phase lose contact in the doping phase diagram toward double- and single- 
layered cuprates. We investigate the mechanism of the disappearance of the coexistence of 
AF and SC in disordered cuprate superconductors in order to understand these doping phase 
diagrams. In single- and double-layered cuprates, electrons on the Cu02 plane experience the 
disorder effect through inhomogeneity in the charge reservoir layer. These impurity potentials 
can be effectively enhanced toward the underdoped region by the effect of many-body correc- 
tions that involve an increase in charge susceptibility. As a result, strong disorder effects are 
expected particularly in the competing regions of AF and SC, where the coexistence phase of 
AF and SC is extremely suppressed. We show the validity of this suppression mechanism by 
considering the Aslamazov-Larkin-type vertex correction to the effective impurity potential 
in the effective mean-field phase diagram. 
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1. Introduction 

The general structure of doping phase diagrams in the underdoped region is one of the 
critical issues with high-temperature cuprate superconductors. In single-layered cuprates such 
as L^-xSr^CuC^ (LSCO), which have been studied extensively since its discovery, it has 
been known that the antiferromagnetic (AF) phase and superconducting (SC) phase do not 
coexist and are separately located in the hole-doping phase diagram while sandwiching a 
spin-glass phase. 1 However, it has recently been discovered by nuclear magnetic (quadrupole) 
resonance experiments that a five-layered cuprate superconductor shows the coexistence of the 
AF and SC states in the underdoped region. 2-5 The discovery of such coexistence indicates 
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that a pure C11O2 plane shows the coexistence of the AF and SC phases, and that the single- 
layered cuprates such as LSCO lose their coexistent phase because of disorder from the charge 
reservoir layer adjacent to the Cu02 plane. The existence of the coexistent phase has been 
substantiated by numerous theoretical approaches on the basis of the Hubbard model 6-14 and 
t-J model. 15-18 However, there has been no established theory for describing the disappearance 
of the coexistence due to inhomogeneity or some other factors on a unified picture. 

The existence of disorder effects from the charge reservoir layer, the so-called out-of- 
plane disorder effects, has been reported both theoretically and experimentally. The impurity 
potential from the Sr site in LSCO is estimated by the local spin-density approximation +U 
(LSDA+U) method. 19 LSDA+U studies have shown that the in-gap band with an oxygen-like 
character shifts further into the gap by ~ leV through a distortion of the Cu02 plane. In 
comparison with the hopping energy estimated as td p ~ leV, this impurity potential through 
the lattice distortion induced by the out-of-plane disorder is not strong enough to reach 
the unitarity limit but is also ineligible for the strong disturbance of the electronic state. 
This situation is different from an in-plane disorder such as the Zn substitution of the Cu 
site where Zn itself becomes a stronger scatterer. In the case of the in-plane disorder, the 
suppression of the SC transition temperature T c in the underdoped region is sharper than 
that in the overdoped region. 20 The suppression of the SC transition temperature T c due to 
an out-of-plane substitution has been observed in single-layered cuprates. 21-24 Its suppression 
of T c occurs with a smaller residual resistivity in the case of the out-of-plane disorder than 
in the case of the in-plane disorder. 21 The pseudogap temperature T* is robust against the 
out-of-plane disorder in contrast to the T c being sensitive to the out-of-plane disorder. 23 

In the double-layered cuprate YF^CusOe+x (YBCO), whose Cu02 plane is less disor- 
dered than the LSCO but not so clean as the five-layered cuprate, the coexistence of a spin 
density wave and superconductivity has recently been observed by neutron scattering. 25 Fig- 
ure 1 shows a schematic phase diagram of the five-layered, double-layered, and single-layered 
cuprates by taking the disorder strength of the Cu02 plane into account. We illustrate in 
their phase diagrams that the AF state, SC state, and their coexistence state (AF+SC) are 
suppressed around the competing regions of the AF, SC, and AF+SC states accompanied by 
the effect of impurity potential. In LSCO, not only does the AF+SC phase vanish but also 
the AF phase and SC phase lose contact with each other because of the stronger impurity 
potential resulting from the substitution of La by Sr. However, the difference between the 
phase diagrams shown in Fig. 1 does not seem to have been clarified yet on the basis of the 
microscopic theory. In this paper, we propose a possible mechanism to describe the difference 
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Fig. 1. Schematic doping-temperature phase diagram for different disorder strengths we propose. 
Antiferromagnetism (AF), superconductivity (SC), their coexistence (AF+SC), and spin-glass 
(SG) phase are drawn. Here, the stripe phase at 1/8- filling and incommensurate SDW are not 
drawn for simplicity. We discuss the difference in phase diagrams among the cuprates from the 
viewpoint of the strength of the disorder potential on the Cu02 plane. 



in doping phase diagrams among clean and disordered cuprate superconductors. 

This paper is organized as follows: In §2, we introduce the concept of renormalized impurity 
potential and consider its behavior in cuprate. In §3, we discuss the role of AL-type correction 
in charge susceptibility and impurity scattering. Then, on the basis of the microscopic theory 
starting from the Hubbard model, we analyze the effects of the Aslamazov-Larkin (AL) term 
on charge susceptibility (§4), and on the phase diagram of disordered cuprates (§5 and §6). 

2. Relationship of Fermi Liquid Corrections between Charge Susceptibility and 
Impurity Potential 

A nonmagnetic impurity potential generally suppresses antiferromagnetism and unconven- 
tional superconductivity. The impurity-concentration dependence of the SC transition tem- 
perature T c is well described by the well-established Abrikosov-Gor'kov formula in many 
superconductors. 26 

In this section, we introduce the concept of the effective impurity potential renormalized 
by electron-electron correlation effects. This concept paves the way for discussing disorder 
effects such as the suppression of AF and SC in cuprates, which are typical strongly correlated 
electron systems. 

The renormalized impurity potential v is defined by considering the vertex correction as 
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Fig. 2. Diagrammatic representation of relationship between the effective impurity potential v and 
the bare potential v through the correction from the vertex function T. 



with 

A(k, q, ie n ) = 1 + T G(k', ie' n )G(k' + q, ie' n )T(k', ie' n , k, ie n , ; q, 0), (2) 
k',n' 

where vq and A are the bare impurity potential and its renormalization factor, respectively. T 
denotes a four-point vertex function. Figure 2 shows the diagrammatic representation of the 
renormalized impurity potential v. The renormalization factor of the impurity scattering for 
the quasi-particle at the Fermi surface is given in the forward scattering limit q = 0: 27 

limA(k F ,q,0) = i T ^, (3) 

where 1/2 = 1 — <9£ R (kp, e)/de\ e= o is the mass renormalization factor and Ffi is the isotopic 
and spin-symmetric component of the Landau parameters. 

A similar enhancement factor can be seen in the charge susceptibility given by the linear 
response theory as 

Xc (q, »w m ) = - dr (n q (r)n_ q (0)) e^ T , (4) 
J o 

where n q (r) = ^ k a cj c+qCT ('r)ci <;o -(T) is the charge density operator, cj^ and Cko- are the cre- 
ation and annihilation operators for the particle with the momentum k and spin a =t,i, 
respectively. We can define the charge compressibility k = dn/dp, which is the uniform limit 
of the static charge susceptibility: 

k= limxc(q,w = 0). (5) 

In a noninteracting system, the charge compressibility reaches kq ~ 2Ny where is the 
density of states at a Fermi level. The relation of the charge compressibility k with k can be 
derived from Fermi liquid theory: 28 

- = (6) 

where m*/m is the mass enhancement factor, which is generally different from 1/z by the 
A:-mass appearing from the k-dependence of the self-energy of electrons. The fe-mass can- 
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eels the mass renormalization factor 1/z. In Lai—^Sr-^MnC^, for example, the effective mass 
m* remains unrenormalized, while the Drude weight D decreases toward the metal-insulator 
transition. At the transition, the Drude weight D = irn*e 2 /m* becomes zero owing to the dis- 
appearance of the effective carrier density n* = n(l + F^/3). 29 However, if the enhancement 
of A and k originate from the anomaly in the Landau parameter Fq — >■ — 1, we can identify 
the origin of the enhancement in the charge compressibility as that in the effective impurity 
potential in the forward scattering limit. 

In the presence of quantum critical fluctuations, vertex corrections affect the enhance- 
ment in the effective impurity potential and residual resistivity. 30-32 As discussed above, the 
enhancement (or suppression) factor of the impurity potential is related to that of the charge 
compressibility. For these reasons, we can expect a strong disorder effect in the material in 
which a large enhancement in the charge compressibility is observed. 

The charge compressibility, or uniform charge susceptibility, has been studied both theo- 
retically and experimentally in cuprate superconductors. Pioneering theoretical works on the 
charge compressibility around the metal-insulator transition have been performed with the 
use of the quantum Monte Carlo method.. 33-36 They have shown in the two-dimensional Hub- 
bard model that the charge compressibility diverges as the doping rate to the Mott insulating 
phase decreases as 5 = 1 — n — > 0. The charge compressibility behaves as k oc 1/5 around the 
Mott insulating phase. This result implies that the transition from a metal to a Mott insulator 
induced by controlling the filling is of the second order. Furukawa et al. have also concluded 
that the divergence k involves the behavior of the Landau parameter, Fq — > — 1 as 5 — > 0. 
Then, they have speculated that the effective mass m* diverges. The 1/5 dependence in k has 
also been obtained from the path-integral renormalization group (PIRG) method. 37 Since the 
quantum Monte Carlo and PIRG methods are applicable to systems with up to about 100 
sites, the appearance of the singular charge susceptibility is not limited in small size systems. 

The experimental study of the charge compressibility has been conducted by photoemis- 
sion spectroscopy. 38-42 The theoretical prediction of the chemical potential shift Afi oc 5 2 , 
which yields k = —dS/dA/u, oc 1/5, is satisfied in LSCO, Bi-2212, and YBCO within the range 
of experimental error. 

Although half-filled cuprates are charge transfer insulators with a different nature from 
the Mott-Hubbard insulator appearing in the Hubbard model, the enhancement in the charge 
susceptibility is observed in some cuprates. Their results suggest that a weak impurity scat- 
tering potential should be enhanced to the extent of the unitarity limit in the underdoped 
region of LSCO, Bi-2212, and YBCO, where an enhanced charge compressibility is observed. 
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Fig. 3. AL term in the charge susceptibility. Solid lines and wavy lines represent the single-particle 
Green function and spin fluctuation propagator, respectively. 

3. Aslamazov-Larkin Term in Charge Susceptibility and Impurity Potential 

In this section, we introduce an approach to describing the singular behavior of the charge 
susceptibility on the basis of the mode-mode coupling theory. 

The divergence of the charge susceptibility in the Hubbard model is known to involve the 
divergence of the AF spin correlation length from quantum Monte Carlo calculations. 34 ' 35 
The singular behavior of the charge susceptibility has also been obtained in a one-dimensional 
Hubbard model by the cellular dynamical mean field theory (DMFT), but this singularity 
cannot be obtained by the single-site DMFT. 43 This suggests the importance of an off-site 
correlation for the appearance of the singular charge susceptibility. There is little doubt that 
the appearance of the singular charge susceptibility in the Hubbard model is closely related 
to the Mott transition. However, we cannot exclude the possibility of the development of AF 
correlation as the origin of the singular charge susceptibility with this background. 

It is desirable to understand the mechanism of the enhancement in the charge susceptibil- 
ity from the structure of the single-particle excitation spectrum in the Green function G and 
the four-point vertex function V in eq. (2). From this viewpoint, it has been known that the 
contribution of the vertex correction shown in Fig. 3, which is the so-called AL term, gives 
a singularity to the charge susceptibility around the AF critical point. 44 The AL-type contri- 
bution to the static charge susceptibility constructed by the two-spin-fluctuation propagator 
X s (q, iw m ) is given by 

AXc,AL(q) = 3T^ [7 3 (q,q / ;^ m )] 2 Xs(q / ,^ m )xs(q / + q,iw m ), (7) 

q',m 

where the factor 3 corresponds to the number of spin directions. The effect of the AL-type 
diagram has been discussed in light of the problem of the electric conductivity anomaly in 
the presence of strong SC fluctuations. 45 73 denotes the coupling strength among one charge 
fluctuation field n(q) and two spin fluctuation fields Si(q) (i = x,y,z), and consist of three 
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Green functions G's for electrons: 

73(q,q';«Wm) = A 2 T^G(k,ie^G(k + q,ie^G(k + q^e n -Hw m ), (8) 

k,n 

where A is the cubic coupling constant in the vertex with the structure of J2k q ^i(l) c {k)c(k + 
<?)• 

Here, we evaluate the effect of critical spin fluctuations that are singular near the incom- 
mensurate AF wave vector Q* in the form 

where rf s is the distance from the magnetic critical point. The spin susceptibility in doped 
cuprate shows its peaks at four incommensurate ordering vectors Q* = 7r/a(l,l±e), n/a(l ± 
e,l). 

Since the dominant contribution to kalCq) arises from the summation at oj m = and 
around q' ~ Q* in the case of rj s <C 2CttT, the contribution of the coherent part denoted by 
eq. (9) is estimated as 

AxcAL(q) = 3T^ 7 3( q ,Q;,0) 7 ,( q ,Q-;0)X: ~ + ^ _ q.,, A(q + q , _ W (10) 

We can see that Ax c ,AL(q) shows its peak at q = Q* — Q*. The maximum is attained at 
q=0: 



Ak al = Axc,al(0) 



= 3x4[ 73 (0,Q*;0)] 2 r^ 



N? 



F 



(2D, T > r] s /2Cir) 

" h t%? )? ■ ffi (3D, T » ' 
In the two-dimensional square lattice with only the nearest-neighbor hopping, 73(0, Q*;0) 
takes a nonzero value except for the half-filling (Appendix A). 

At low temperatures, % » 2CirT, Akal shows a less singular behavior because the finite- 
u m contribution is not negligible. At zero temperature T = 0, the frequency summation in eq. 
(7) can be approximated by the integration with respect to uj. As in the case of % » 2CttT, 
the maximum is attained at q = 0: 

roc dw iV| 



r 

Aacal ^ 3x4[ 73 (0,Q*;0)] 2 -2V / 

~7 Jo 



2tt fa + Aq' 2 + Coj) 2 
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- 24M0,Q-;0)]sfE^ 



q 



— < 



^Mlog^ (2D,T = 0) 

12 [7 3(0,Q*;0)]^^_^^ (3Dr = 0) ' (12) 



7T 2 AC 

where q c ~ 1/a is the short-wavelength cutoff and of the order of the inverse of the lattice 
constant a. The singularity in Akal at quantum critical points is weakened in two dimensions 
or suppressed to the cusp singularity in three dimensions. The charge susceptibility shows the 
power-law singularity near the magnetic instability 7] s ~ at finite temperatures, where the 
critical exponent of Akal are the same as that of the SDW susceptibility Xs(Q*>0) in two 
dimensions, and half of that of Xs(Q*>0) in three dimensions. However, this critical exponent 
is not always reliable because it is easily changed by other processes such as the mode-mode 
coupling effect among charge-density channels. 

Other peaks in Ax c ,AL(q) appear at q = Q* — Q* / 0. These peaks are generally smaller 
than Axc,al(0) because 73 is smaller for such mode coupling processes. The case of Q* = — Q* 
corresponds to a stripe charge fluctuation. The stripe charge fluctuations are known to be 
widely observed in cuprate superconductors. Moreover, a stripe ordered phase is stabilized 
instead of the SC phase at a doping rate 5 = 1/8. We note that Ax c ,AL(q) shows the emergence 
of stripe fluctuations, but unfortunately, is not accountable for the stabilization of the stripe 
ground state with completely suppressed superconductivity only at S = 1/8 in cuprates. 

Here, we consider the effect of the AL process on the impurity scattering problem. By 
taking into account the AL-type correction to the impurity scattering process, we can obtain a 
singularly strong disorder effect around AF critical points. Such a correction must be expressed 
in diagrammatic form as shown in Fig. 4. 

When the charge susceptibility is nearly critical at q = 0, the renormalized impurity 



8/32 



J. Phys. Soc. Jpn. Full Paper 

potential v at the Fermi level can be expanded around the forward scattering as 

s(q) ~ ^rkF'°° (q) ' (13) 

where tj c , the distance from the critical point with respect to a uniform charge density fluctua- 
tion, corresponds to the square inverse of the charge density correlation length r/ c oc This 
kind of vertex correction leads to a long-ranged impurity potential even if the bare potential 
vo is short-ranged. The range of the impurity potential is given by £ c , which diverges at the 
critical point. 

We discuss the scattering rate of quasi-particles and the residual resistivity due to renor- 
malized impurity scattering under the condition of a spherical Fermi surface. The scattering 
rate 1/r, which causes the suppression of the superconductivity and antiferromagnetism, re- 
lates to the imaginary part of the self-energy. It is attained for the impurity concentration 
C; mp within the Born approximation as 

- ~ 27r Cimp ^|i)(k-k / )| 2 z<5(r(k / )) 
Tk k' 

= 2vriV FCimp (- \ , (14) 

where £*(k) is the quasi-particle dispersion and is the angle between k and k'. (■•■)ps 
denotes the average on the Fermi surface. Here, note that the Born approximation is not 
applicable in the strong impurity scattering region around the criticality rj c ~ 0. In such 
cases, eq. (14) should be replaced with 1/rk = 2ci mp Imi (k, k) with the use of the t-matrix 
t(k, k') = v(k— k') + ^ k » v(k — k")GR,(k", 0)t(k", k'). If we ignore the momentum dependence 
in v(k — k'), the t-matrix is approximated as i _1 = v^ 1 — iirNp. By this reason, the scattering 
amplitude saturates to 1/r ~ 2ci mp /irNp toward the unitarity limit. However, we can examine 
the feature of the enhancement even if the Born approximation is adopted. By taking the 
average in eq. (14), the behavior of the scattering rate near the criticality rj c ~ is 

l/r?c /2 (2D) 

• (15) 

1/ Vc (3D) 



1 

oc < 

Tk F 



The electric resistivity shows a different behavior from the scattering rate. The residual 
resistivity po due to the impurity scattering is obtained as 



1 

po oc — 

Ttr 



Vq(1 — cos ( 
(r ]c + 2Akl(l-cos0)y , ,, s 
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The contribution of the forward scattering is subtracted in po so that the lifetime in transport 
phenomena, r tr , is different from in general. We can easily take this average. As a result, 
we obtain 

Vv^fc (2D) 



Po oc < 



(17) 

log(l/»fc) (3D) 



The residual resistivity also shows a divergence at critical point rj c = 0. However, its power 
of divergence is weak compared with that of the scattering rate 1/t^ f . Thus, even if the 
enhanced impurity potential, due to critical charge fluctuations, destroys the long-range order, 
the residual resistivity does not show a visible enhancement in certain cases. Therefore, even 
though the residual resistivity po is insensitive to the impurity potential in contrast with the 
quasi-particle lifetime r^, the SC transition temperature T c can vanish at a smaller residual 
resistivity in the underdoped region than in the overdoped region. 

4. Aslamazov-Larkin-type correction to the charge susceptibility in the Hubbard 
model 

In this section, we analyze the role of the AL-type vertex correction in the Hubbard model 
with the use of the spin-fluctuation propagator at the level of random phase approximation 
(RPA). The consideration of the spin susceptibility within the RPA level is consistent with 
the analysis of the mean-field phase diagram for the antiferromagnetism. For comparison, we 
also study the effect of the Maki-Thompson (MT) term consisting of aq diagram simpler than 
the AL term. 

The two-dimensional Hubbard model in a square lattice is described by 

H = Yj £ ( k ) C L C k<x + U Yl C k+qf C kt C k'+ C k'+q+ • (18) 
k,<r q,k,k' 

The dispersion in cuprates is well fitted by e(k) = — 2t(cos k x + cos k y ) + At 1 cos k x cos k y with 
the hopping integral between the nearest-neighbor sites t and the next-nearest neighbor sites 
t'. In this study, we choose t' = for simplicity and owing to the limitation of RPA analysis 
(see §6). In this section, we perform a diagrammatic analysis by putting the Green function 
for noninteracting electrons G^\\s.,ie n ) = [ie n — e(k) + p]~ l into all G a (k,ie n ). This means 
that no self-energy correction is considered in the Green function G. 

The charge susceptibility is expressed in terms of the vertex function F and the Green 
function G as shown in Fig. 5, and its analytic expression is given as 

Xc(q) = 2 X0 (q) + T 2 £ G a (k)G a (k + q)T aa i (k, k'; q)G a i (k')G a i (k! + q), (19) 

k,k',cr,cr' 
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Fig. 5. Diagrammatic representation of relationship between charge susceptibility and vertex func- 
tion. 



(a) 



(b) 



- a 



Q 

o 
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Fig. 6. Maki- Thompson-type vertex correction in RPA. Broken lines show the on-site interaction U. 
Each graph shows the contribution of (a) the transverse spin fluctuations and (b) the longitudinal 
spin fluctuations and charge fluctuations. 



where we have introduced the abbreviations k = (k, ie n ) and q = (q, iuj m ) with the Matsubara 
frequencies, e n = (2n+ 1)ttT for electrons, and u m = 2irmT for bosonic excitations. We define 
the irreducible particle- hole bubble xo{q) as 

X0 (q) = -Tj2G(k)G(k + q), (20) 
k 

where the frequency summation is analytically performed in the non-interacting system with- 
out disorder as 

Xo(9) = E- fj 7VpYZ ( 21 ) 

^ iu m - e(k + q) + e(k) 

where /k = [exp((e(k) — n)/T) + is the Fermi distributi on function. 

We evaluate the particle-hole irreducible vertex function T\ representing MT processes and 
AL processes using spin fluctuations within RPA. The two MT terms in the vertex function 
are shown in Fig. 6: 

iYmtw = E r w> ( 22 ) 

i=a,b 

tV-w* - (23) 
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Fig. 7. AL-type vertex correction in random phase approximation. Even (odd) denotes the number 
of particle-hole bubbles xo- 



Six AL terms shown in Fig. 7 are given as follows: 
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The second-order diagram in U has been subtracted from Ty -C '"' ' for the purpose of avoiding 
double counting of the MT term i\mt- The AL terms consist of the particle- hole bubbles for 
j-i(»-c,-,/)^ w hj c h i s the contribution of the longitudinal spin fluctuations and appears when 
two propagating fluctuations have the same parity with respect to the number of particle-hole 
bubbles, and the particle-hole ladder for Y^~ g ' h \ which represents the contributions from the 
transverse spin fluctuations. Their total irreducible vertex function is expressed as 

^ l,aa 

(k,k';q) = US^+Ti ,MTcrcr' (k,k';q)+T 1 (k,k';q). (32) 

We numerically investigate the role of the MT and AL terms for the charge susceptibility. 
To this end, we calculate the total charge susceptibility using the relation 

XM = 2xo(9) - 2U X 2 o(q) + AxMT(g) + AxaUq), (33) 

where the third and fourth terms are respectively defined as 

A X mt(</) = T 2 G r7 (k)G a (k + q)r 1MTaa/ (k,k';q)G a/ (k')G af (k' + q), (34) 

k,k',cr,cr' 

Axal(^) = T 2 G c {k)G a {k + q)T lAlj(TC ,{k,k'-q)G a ,{k')G c ,{k' +q). (35) 

k,k' ,cr,cr' 

For the computation of the MT- and AL-type contributions to the susceptibility, the fast 
Fourier transformation (FFT) algorithm is very useful from the viewpoint of numerical effi- 
ciency. The application method of FFT to the AL term is summarized in Appendix B. We 
perform the calculations under the condition of 32 x 32 lattice sites with 512 Matsubara 
frequencies. 

Figure 8 shows the q dependence of each term in the static limit (w m =0). There are several 
peaks in the AL term AxAL(q), as discussed in §3. The maximum is attained at q = 0. The 
positions of the peaks in the MT term A%mt almost coincide with the incommensurate 
ordering vector Q*, peak positions of Xo(q)- Figure 9 shows the doping dependence of each 
term. Axal rapidly develops toward the AF transition, as predicted in §3. In contrast, Axmt 
gives only a small contribution to Xc even at q 7^ 0. 

The singularity in the AL term Axal, discussed so far, arises through the divergence of 
m* in eq. (6). Hereafter, we show that the repetition of the AL-type vertex correction can 
give the enhancement of the factor (1 + Fq) -1 . The actual calculation for solving the Bethe- 
Salpeter equation based on the irreducible vertex is generally difficult for arbitrary diagrams, 
but we can evaluate the summation of several terms in the lowest-order terms. To check the 
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2X0 (q) 




Fig. 8. q dependence of the Aslamazov Larkin term Axal(q), the MT term Axmt(q), and the 
noninteracting part 2xo(l) m the static charge susceptibility for a parameter set U = 2t,n — 0.875, 
T = 0.05*, and 1 - U X a(Q*, 0) = 0.015. 
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Fig. 9. (Color on line) Doping dependence of the AL term (red), maximum of MT term (green), and 
non-interacting part (blue) in the charge susceptibility when U = 2t and T = 0.05*. AF denotes 
the AF-ordered region. 
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Fig. 10. (Color on line) Higher-order effects on AL term (upper panel) and MT term (lower panel) 
for the charge susceptibility at q = 0. Both terms are evaluated for the parameter set U = 0.05i 
and T = 0.05i. 



higher order effects in the MT and AL terms, we consider the following correction in the 
charge susceptibility summed up to the n-th order of the irreducible vertex: 



AxMT.n = GGTi mtGG H + GG (T^utGG) , 

AxAL.n = GGT 1AL GG + ■ ■ ■ + GG (T hAL GG) n . 



(36) 
(37) 



We do not show the higher-order terms including both I\mt and I\al for simplicity. 

Let us see the tendency of the behavior in the n — > oo limit of A^AL.n- Then, Akal m 
eq. (11) and eq. (12) in two dimensions are replaced with a geometric series: 44 



Akal -> Axal,oo(0) 



Ak A l/ 
Ak A l/ 



1 | 74 N*T 



1 + 



. log ^ 



(T > rj s /2CTr) 
(T = 0) 



(38) 



it AC ""° r/s 

where 74 is the quartic mode- mode coupling constant of spin fluctuations. 74 for the quartic 
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coupling among four fluctuation modes near the q ~ Q* is expressed as 

74 = A 4 T^G 2 (p,*e n )G 2 (p + QVe n ). (39) 

p,n 

74 takes a negative value at T < 0.524|/i| in the case of a two-dimensional square lattice (see 
Appendix C). At that temperature, 5kal is enhanced in a region far from the AF critical 
point. 

Figure 10 shows the behaviors of Axmt> and Axal,« for n = 1, 2, and 3. The singularity 
in Axal(0) increases as the order n of the diagrams increases. These results show that AL- 
type vertex corrections give a negative value to the Landau parameter Fq as well as 74 < 
around the AF criticality in the low-temperature region. Concerning the MT term, A^mt does 
not show a remarkable change with higher-order terms. The higher-order terms in the MT 
terms are negligible, but those in the AL terms considerably increase the charge susceptibility 
around the AF criticality. 

Note that, throughout this section, we have evaluated the anomaly in the charge compress- 
ibility with the use of the Green function in a noninteracting system. Namely, the self-energy 
correction arising from electron-electron interaction in the Green function is not considered. 
To improve this point, it is desirable to take self-energy correction consistent with vertex 
correction into account so as to satisfy the conservation law. 46 However, conserving approxi- 
mation such as fluctuation exchange approximation 47 do not correctly produce a coherent part 
in the single particle Green function and is also not capable of describing the Mott transition. 
Since the singular charge compressibility originates from the formation of the gap in single 
particle excitation at half-filling, self-energy correction has to reproduce the appearance of 
the charge gap. However, it is still challenging to construct such a conserving theory. For this 
difficulty, we ignore self-energy correction for simplicity and consider only vertex correction 
to the impurity potential in the next section. 

5. Disappearance of coexistence of antiferromagnetism and superconductivity 
due to large impurity potential enhanced by AF fluctuations 

In the remaining sections, it is shown that an inhomogeneity leads to the suppression of 
both the Neel temperature Tn and the SC transition temperature T c . To this end, we extend 
the theory of the charge susceptibility taking the effective mass enhancement due to the AL 
term to the impurity scattering problem in cuprates. As discussed in previous sections, we 
focus on the problem of how the impurity scattering potential enhanced by the singular charge 
fluctuation affects Tn and T c . Here, we introduce the two-dimensional Hubbard model with 
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Fig. 11. Self-energy diagram within Born approximation with respect to renormalized impurity po- 
tential. The shaded triangle represents the renormalized impurity potential. 

the d x 2_ y 2-wave pairing interaction and the impurity potential vf. 

H = e ( k ) c L c ka + u ^2 n it n ii 

kcr % 

(40) 

k ia 

where the form of the gap function 0(k) = cos k x — cos k y is chosen to be consistent with 
observations in cuprates. The bare impurity potential is given by Vi = vq for the impurity 
sites, and Vi = for the other sites. For L^-^Sr^CuO,!, it corresponds to 

vq at Cu02 sites adjacent to Sr, 

(41) 

at Cu02 sites adjacent to La. 
Namely, we consider a situation where the impurity concentration and the doping rate 5 take 
the same value. Of course, there exist other types of inhomogeneities, such as the inhomo- 
geneity of hopping integrals and the inhomogeneity of interactions on the Cu02 plane, even 
if they are caused by the out-of-plane disorder. However, we assume here the local impurity 
potential in order to study the qualitative aspect of impurity potential renormalization. 

Hereafter, we discuss the effects of vertex correction due to AF critical fluctuations neglect- 
ing vertex corrections from the d-wave attraction U^. We define the renormalized impurity 
potential v within the first order approximation in the irreducible vertex T\ as 

6(k,k + q,ie n ) = v A(k, q, ie„), (42) 

with 

A(k, q, ie n ) = 1 + T G(k', ie' n )G(V + q, ie'^T^k', ie' n , k, ie n , ; q, 0). (43) 

k',n' 

Hereafter, we use the irreducible vertex function Ti given by eq. (32), which includes magnetic 
critical fluctuations within the RPA. 

We evaluate the self-energy in the single-particle Green function within the self-consistent 
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Born approximation. With the use of the renormalized impurity potential, we obtain the 
self-energy shown in Fig. 11 as 

S imp (k,ze n ) = c imp J^V(k,k',ie n )G(k',ie n ), (44) 
k' 

where c; mp is the impurity concentration c- nnp = |1 — n\. Here, we have assumed that the bare 
impurity potential is a real number so that the renormalized potential v satisfies v(k, k', ie n ) = 
v(\s!, k, ie n ). Even if the bare potential is real, the renormalized impurity potential can be a 
complex number in general. Even so, the imaginary part of the analytically continued self- 
energy still gives the damping rate of the quasi-particle as l/2rk = — Im£i mp (k, e + iS). 
We must comment on the first-order term of v in the self-energy, which is given by 

A£« (k, ie n ) = c imp £(k, k, ie n ). (45) 

For the bare potential, since it has the general form of t>o(k, k', ie) = t>o(k — k'), the first-order 
term is constant, giving only a chemical potential shift. Unlike the bare potential, eq. (45) 
gives the k and e n dependences. Unfortunately, the sign of the renormalized potential v in Eq. 
(45) depends on the sign of the bare potential vq. Thus, it shows an anomalous behavior that 
the damping rate easily becomes negative. However, this term must give only the chemical 
potential shift in the Green function G constituting the self-energy like 

£(k,ie n ) = Tj]F(q,iw m )G(k + q,i€ n + ia; m ), (46) 

q,m 

V(q,iu m ) = U+-- —— : r + - - : - U z x Q {^iu m ). (47) 

2 1 - Uxo(<l,iu m ) 2l + Uxo(q,iu m ) 

This self-energy S is related to the irreducible vertex function by = ^i(k, k'\ q = 0). We 

do not consider this kind of self-energy correction in the present paper. However, because eq. 

(45) should originally be absorbed into eq. (46), we eliminate this first-order term. 

We note that the Born approximation is valid in the case of a dilute concentration of 
impurities with a weak potential. In addition, since it overestimates the damping rate for a 
strong impurity potential, the application of the t-matrix theory is desirable in such a case. 
Nevertheless, here we use the Born approximation for simplicity because better theories, such 
as the t-matrix approximation, are numerically difficult to adopt for renormalized impurity 
potential including complex k and k' dependences. 

If we have the calculated self-energy Si mp (k, ie n ) due to impurities in eq. (44), the single- 
particle Green function is obtained from the Dyson equation: 

G _1 (k, ie n ) = ie n - e(k) + fi - S imp (k, ie n ). (48) 
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Using this type of Green function, we repeatedly calculate the vertex function in eq. (32), the 
renormalized impurity potential in eq. (42), and the self-energy in eq. (44). This self-consistent 
calculation scheme is performed until G, Xi mp and T converge. If a self-consistent solution is 
obtained, we can evaluate the effect of renormalized impurity scattering on the SC transition 
temperature T c and the Neel temperature Tn. Here, the carrier density n is evaluated as 

n = 2T^G(k,ie n )e ie " <5 

k,n 

= 1 + 4T ReG(k,ie n ). (49) 

k,n>0 

We note here that dn/djJL calculated from the carrier density in eq. (49) and from the Green 
function in eq. (48) does not give divergence toward AF criticality because G(k,ie n ) does 
not include the self-energy correction of eq. (46) with an effect of critical AF fluctuations. In 
this sense, relation eq. (5) is broken due to difference in the level of perturbation between k 
calculated from the /j, dependences of n and Xc constructed by using an approximate vertex 
function T±, eq. (32). 

The AF transition temperature Tn is determined with the use of the converged Green 
function by 

I- a = l/max[xo(Q,0;T = T N )], (50) 
Xo(Q,0;T) = -T^G(k,ie n )G(k + Q,ie n ). (51) 

k,ra 

The ordering vector Q, which gives the maximum Xo(Q>0), is not restricted to the commen- 
surate AF vector Q = (tt, tt) in general. Thus, hereafter, our criterion furthermore gives an 
incommensurate spin density wave transition. The true solution should be given at a = 0. 
However, this criterion never gives the AF and SDW transition because the AL term gives 
the relation rj c = rj s in eq. (14) and the damping rate diverges at the true AF critical point 
rj s = a = 0. This serious problem is solved by the extension tothe f-matrix theory in which the 
damping rate saturates at a constant value at r] c = 0, as discussed in §3. Because the exten- 
sion to the ^-matrix theory of the renormalized impurity potential is very difficult to compute 
within reasonable computation time, we choose the condition of a = 0.005 to estimate Tn- 
This a corresponds to the situation where the 2D AF coherence length takes £af ~ 13a, with 
a being the lattice constant, at U = 2 and in the half-filling. 

The SC transition temperature T c is given by the Thouless criterion 

l = U d T£|G(k,ie n )|V(k). (52) 

k,n 
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Fig. 12. Mesh in q-space for 32x32 lattice for numerical calculations. The q-summation in the Born 
self-energy, eq. (53), is performed over the center in the square surrounded by points with X 
symbols. 

Since the impurity effect is included in the dressed Green function, a solution of eq. (52) gives 
a decrease in T c in disordered systems. 

The numerical summation in eq. (44) requires a much longer time for calculations than 
that for xoi eq. (51). Then, we rewrite Eq. (44) as 

S imp (k, ie n ) = CimpVg ^2 G(k, ie n ) + c imp ^ (v 2 (k, k + q; ie n ) - Vq) G(k + q, ie n ). (53) 

k q 

The first term corresponds to the Born approximation for the bare impurity potential vq, 
and its summation can be easily performed. In contrast, the second term is numerically more 
difficult to evaluate. Here, since the structure of the vertex correction appears at q = 
and q = (vr,7r), we approximate the summation using q-points with a high-priority around 
q = and (tt,it). The resulting mesh points are shown in Fig. 12 for the 32 x 32 lattice. By 
discretizing the q-integral in eq. (53) with the use of the trapezoidal formula in each squared 
region, we can derive the weight of each mesh point in the discretized summation. 

Although we take into account only the self-energy owing to impurities in the present pa- 
per, the above framework enables us to discuss how the mean-field phase diagram is influenced 
by renormalized impurity scattering. 

6. Doping Phase Diagram in High-T c Cuprates with Disorder 

First, let us investigate the behavior of the enhancement factor of the impurity potential 
and the Born self-energy in the limit of v o — > for a parameter set U = 2t and n = 0.875, 
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Fig. 13. Real and imaginary parts of the renormalization factor A(k, k, iirT) for the impurity potential 
as a function of the momentum k of scattering electron. The momentum transfer q is chosen to 
be zero. 



where % = 1 — f7max[xo(Q, 0)] = 0.015. 

Figure 13 shows the incident momentum (k) dependence of the renormalization factor of 
the impurity potential A(k, q, ie n ) in the forward scattering limit q = 0. Its real part exhibits a 
peak at the Fermi surface, but the peak in the imaginary part deviates from the Fermi surface. 
Thus, only the enhanced real part contributes to low-energy impurity scattering processes. 
Note that the renormalization factor in our theory satisfies the relation lim q ^.o A(k, q, + 
iS) = 1 + <9X(k, + iS)/dfj,, where the self-energy due to the electron correlation E(k, ie n ) 
is given by eq. (46). In a normal Fermi liquid, ImS(k, e + id) oc max(e 2 , n 2 T 2 ) is satisfied. 48 
Thus, lim q _;.o ImA(k, q, + i5) = (9ImS(k, + i5)/dfi generally takes a nonzero value at 
finite temperatures. Therefore, the renormalized impurity potential generally takes a complex 
value even in the static limit (ie n — > + i5), even if the bare potential is real. For these 
reasons, the relation v(k,k',ie) = v*(k',k, ie) required in the bare potential is not satisfied. 
Figure 14 shows a momentum transfer (q) dependence of the renormalization factor. Its 
global maximum appears at q = 0, as in the case of the AL-type correction of the charge 
susceptibility. Secondary peaks are also obtained near q = as the AL term in the charge 
susceptibility. Correspondingly, we show the structure of the integral of the self-energy in 
eq. (44). The renormalized impurity potential is mainly enhanced in the forward scattering 
channel. Thus, the Born self-energy £i mp can be approximated near the criticality as 

S im p(k, ie n ) oc v 2 (k,k,ie n )G(k,ie n ) 

~ Re[v 2 (k,k,ie n )]lmG(k,ie n ), (54) 
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Fig. 14. Real and imaginary parts of impurity potential renormalization factor A(k, k + q, inT) as 
functions of transfer q. The incident momentum is chosen to be k = (-7T, n), around which the real 
part shows its peak. 




(0,0) (0,vr) (tt.tt) (0,0) 

Fig. 15. Comparison of real part of v 2 (k, k, inT) (solid line) and ImG(k, inT) (dotted line) in the 
k-space at n = 0.875 and T = 0.05t. The variable k sweeps along the broken line with the arrow 
shown in the inset. The Fermi surface for this case is drawn with a dotted line in the inset. 

where we have ignored the contribution of ReG(k, ie n ) in the second line, because ReG(k, + 
i8) approaches zero at the Fermi surface k = kp. Figure 15 shows the k dependence of 
Re[i) 2 (k, k, iiT)] and ImG(k, iirT). The peak in ImG(k, inT) corresponds to the Fermi surface 
momentum. The square of the renormalized impurity potential v 2 takes both positive and 
negative values. However, since it takes a positive value at the Fermi surface in the (0, ir) 
direction, the scattering rate of quasi-particle l/2rk ~ ImSi mp (k, iirT) is largely enhanced at 
the Fermi surface. 
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Fig. 16. Real and imaginary parts of the enhancement factor of the Born self-energy £i mp (k, iirT) /vq, 
with £i mp given by Eq. (53), in the limit vq — > 0. 

Figure 16 shows the real and imaginary parts of the Born self-energy, eq. (53), in the first- 
order correction with respect to Vq. On the other hand, the Born self-energy without vertex 
correction at n = 0.875 becomes almost independent of k, in the low Matsubara frequency 
limit as 



This has essentially no k dependence. One can see that the scattering amplitude 
ImSi mp (k, iirT), is enhanced around the Fermi surface particularly in the (0,7r) antinode 
direction. Not only the imaginary part but also the real part of ImSi mp (k, iirT) exhibits a 
visible k dependence. Since Ci m pUg is definitely positive, real-part correction works by reduc- 
ing the next-nearest-neighbor hopping t' . The above results display the enhancement of the 
scattering rate due to critical AF fluctuations. 

Let us see the characteristic of the resulting doping phase diagram obtained by taking 
into account the impurity potential renormalized by the vertex correction. Figure 17 shows 
the resulting doping phase diagram in the extended Hubbard model with disorder, eq. (40). 
As we expected from previous discussions, we find a remarkable suppression effect on the 
transition temperatures Tn and T c around the competing region of antiferromagnetism and 
superconductivity. This is due to the enhancement of the impurity potential renormalized 
by AL-type vertex correction. In fact, such a strong suppression of Tn and T c cannot be 
understood without the vertex correction shown in the inset of Fig. 17, for which no such 





(55) 
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Fig. 17. (Color) Doping phase diagram of disordered extended Hubbard model, eq. (40), with vertex 
correction for impurity potential, vo denotes the strength of the bare impurity potential. Each line 
represents the transition temperature of the antiferromagnetism (AF) in red, the incommensurate 
spin density wave (ISDW) in blue, and the superconductivity (SC) in green. Inset: phase diagram 
without vertex correction for impurity potential. 



suppression appears 

In the phase diagram in Fig. 17, the AF phase shows a reentrant feature as the temperature 
is lowered. However, because the coefficient of the quartic term 74 in the Ginzburg-Landau 
free energy takes a negative value in the low-temperature region T < 0.524|/i|, the AF transi- 
tion is expected to be of the first order in the low-temperature region where dT^(n) / 'dn < 0. 
However, in order to study the reality of the reentrant feature or the first-order transition, 
a more reliable theory such as the t-matrix theory or a theory beyond mean-field analysis is 
required. For an incommensurate spin density wave transition, it is difficult to conclude the 
order of the transition because it requires more detailed analysis of free-energy structures con- 
sidering couplings among multiple ordering vectors. The Vq dependence of the doping phase 
diagram in Fig. 17 reproduces the tendency that both antiferromagnetism and superconduc- 
tivity disappear around the competing region of the AF and SC phases when we move from 
five-layered compounds 2-5 to the double-layered compound YBCO. 25 Moreover, the result- 
ing phase diagram reproduces the emergence of the incommensurate SDW phase observed in 
YBCO. 25 However, we must be careful with the application of our theory to more disordered 
cuprates such as LSCO. In the underdoped LSCO, which experiences a stronger disorder, 
the possibility of a spin-glass state, 49 a charge-glass state, 50 and the Anderson localization 
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Fig. 18. (a) Self-energy correction and (b) vertex correction for the spin susceptibility by the renor- 
malized impurity potential. In this study, we ignore the latter vertex diagram (b). 

due to the two-dimensional character should be discussed. Because the Born approximation 
is not suitable for a strongly disordered system, our theory cannot describe such insulating 
characteristics originating from a strong inhomogeneity. Thus, our theory should be applied 
to systems ] having a moderate (out-of-plane) disorder. 

In contrast with the experimental phase diagram, the AF phase is more robust against 
the impurity potential than the SC phase in our results. The actual critical doping rate for 
the disappearance of antiferromagnetism is obtained as S c ~ 0.12 for five-layered cuprates, 
0.05 for YBCO, and 0.02 for LSCO. In our analysis, however, it is obtained as S c ~ 0.07 for 
vq = t where the SC transition temperature T c is lowered below half of the T c for vq = 0. 
There are two reasons for this discrepancy. First, we have ignored vertex corrections for the AF 
susceptibility as shown in Fig. 18(b). Different from that for the anisotropic superconductivity, 
the vertex correction for the AF susceptibility has a role of destabilizing the Neel order together 
with self-energy corrections, as shown in Fig. 18(a). However, we had to ignore this kind of 
vertex correction because of the long computational time required to evaluate the so-called 
diffuson diagram in which vertex corrections are summed up to the infinite ladder for the 
renormalized impurity potential. Second, both the effects of strong spin fluctuations and 
strong pairing fluctuations are missing in our mean-field analysis. Note that cuprates have a 
larger Coulomb repulsion such as U ~ 6t or 8t. In this region, the pairing interaction goes to 
a relatively strong coupling regime. In the presence of atrong pairing interaction, a Cooper 
pair with a short coherence length is realized. In such a case, the SC state becomes robust 
against impurities in general. 51 

Although electrons in cuprates experience a strong Coulomb repulsion, we have analyzed 
in the weak-coupling situation that U = 2t because the mean-field approximation gives a 
greatly overestimated Tn for large-C regions. Our weak-coupling assumption can be justi- 
fied to some extent by regarding U as an effective low-energy vertex function in the spin 
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channel. Even if the bare U is as large as the bandwidth 8t, the lo-energy effective vertex 
function is suppressed to ~ 2t at half-filling in the two-particle self-consistent theory 52 and 
self-consistent fluctuation theory. 53 These theories are attempts to determine a constant U sp 
which is an approximate spin-channel vertex function in a reliable way. The suppression of the 
effective interaction in the spin channel has also been observed in the pseudo-potential par- 
quet approximations. 54 These saturating phenomena of magnetic interaction are well known 
as the Kanamori-Briickner screening effect in the problem of the ferromagnetic transition 
in itinerant electron systems. Owing to the suppression of the effective interaction trigger- 
ing antiferromagnetism, our mean-field analysis is expected to give a qualitatively reasonable 
description for cuprates. 

In addition, we have not considered the next-nearest neighbor hopping t' and the third- 
nearest neighbor hopping t" , which are necessary for reproducing the quasi-particle dispersions 
in cuprates. In weak-coupling cases, however, the effect of t' makes the AF phase appear only 
away from the half-filling. 55 Then, we would have to require a reliable calculation method in 
the large-C/ regime for the purpose of treating an actual band dispersion in cuprates. 

7. Conclusions 

We have investigated the doping phase diagram of high-T c cuprate superconductors with 
a moderate out-of-plane disorder in order to understand the mechanism by which the co- 
existence of antiferromagnetism and superconductivity exists in five-layered cuprates, but is 
suppressed toward double- and single-layered cuprates as the number of layers decreases. We 
hypothesized, in solving this problem, that the impurity potential is enhanced by the many- 
body corrections in the underdoped region, as shown in Fig. 1. In fact, the experiments on 
the charge susceptibility and discussion on the basis of the Fermi liquid theory suggest that 
the enhancement of the impurity potential should occur in the underdoped cuprates through 
the many-body effects leading to the enhancement of the charge susceptibility. 

On the basis of the microscopic theory taking the scattering process through the renormal- 
ized impurity potential into account, we have shown that the effective impurity potential can 
actually be enhanced in the underdoped region by the AL-type process which also causes the 
enhancement of the uniform charge susceptibility. Although the AL-type corrections to the 
impurity potential gives a 1/z factor, they lead to the remarkable effects that the coexistence 
region of AF and SC decreases in size and that the peak of the SC transition temperature 
T c moves away from the AF critical point with increasing the disorder strength in the doping 
phase diagram. 
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As can be seen in Fig. 17, the enhancement of the disorder effect is actually substantiated 
within the analysis by the mean-field approximation and Born approximation. The resulting 
doping phase diagram can account for the phase diagram of YBCO, but not for the emergence 
of the glassy state in LSCO, relatively strongly disordered cuprates. This problem stems from 
the limited available range of the Born approximation, however, it allows the validity of the 
concept of the renormalization of disorder potential due to strong electron correlations. 
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Appendix A: Cubic mode-mode coupling strength among charge fluctuations 
and AF fluctuations 

The coupling strength among one charge fluctuation field at q = and two AF fluctuation 
fields at Q* = (ir, it) is given by 

73 (0, Q* ; 0) = T (-■ 7^—) 2 -■ nTTTvTX • ( A4 ) 

^ \te n - e(k) + nj ie n - e(k + Q ) + fi 

Here, we suppose that the single-particle dispersion in the square lattice is given as = 
— 2t(cos k x + COS k y ). This dispersion satisfies the perfect nesting condition £k+Q = — e^. Then, 
the coupling 73 is rewritten as 

/ l \2 x 



73 ~ T y Np / de 



00 



ie n - e + fi J ie n + e + fj, 



iirN F T 



2 

n>0 



(ie n -fif {ie n + ii) 2 _ 



= -2kN f Th^ —^ (A-2) 

n > ( £ n + ^) 

In the region \fi\ <C ttT, 73 can be expanded by //. The leading term is obtained as 

2N F fi ^ 1 



73 



7T 2 T 2 ^ (2n + l) 3 

n>0 v ' 
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which is linear in fi. In contrast, in the region 3> irT, the summation can be approximated 
by an integral with respect to e as 



f°° e 
73 ~ —Npfi / de 

Jo 



(e 2 + /i 2 ) 2 

= (A-4) 

The right-hand side of eq. (A-2) is equal to (iV F /47rT)InT0'(l/2 + i^/2vrT), where V' is the 
tri-gamma function. Therefore, I73I reaches its maximum at ji ~ /U max , which corresponds to 
the condition ReV'*- 2 * 1 (1/2 + i[i max /2irT) = 0, ?/>( 2 ) being the tetra-gamma function, because 
the following identity holds 



Then, |/x max | is given by |// ma x| — 1-91T, (See Appendix C). 

Appendix B: Numerical procedure for calculating Aslamazov-Larkin term 

The FFT algorithm is useful in the evaluation of the AL term where the convolution 
integrals are included. In this appendix, we explain how it is applied to the calculation of 
AL terms. We consider the following AL term in the renormalization factor for the impurity 
potential: 

A A L(k,ie„;q) = T 2 ^ G(k', <)G(k' + q, <)G(k' + q', ie' n + iu m ) 

k',ie^,q',iw m 

xr s (q, iu m )T s {<\ - q, iuj m )G(k + q, ie n + iw m ). (B-l) 

It is more convenient to directly calculate the three-point vertex function Aal than the four- 
point vertex T^al because it is expressed by fewer momentum indices. The integral in the 
AL term can be rapidly performed with the use of the FFT at a fixed q. For this purpose, we 
disassemble eq. (B-l) as 

A(k,ie;q) = G(k, ie n )G(k + q, ie n ), (B-2) 

B(q',io; m ;q) = ^ A(k, ie n ; q)G(k + q', ie n + iu m ), (B-3) 

C(q',zu; m ;q) = B(q, iu: m ; q)r(q, iuj m )T(q - q', iu m ), (B-4) 
A A L(k,ie n ,q) = ^ C(q', iw m ; q)G(k + q', ie n + iw m ). (B-5) 

We can obtain the renormalization factor Aal through eqs. (B-2) - (B-5). Note that eqs. (B-3) 
and (B-5) are expressed as the convolution form. Thus, its numerical calculation requires only 
the number of times of 0(N log N) for a fixed transfer momentum q. Moreover, we can obtain 
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the AL-type correction for the static charge susceptibility from the three-point function Aal 
as 

Axc,AL(q) = ^2 A A L(k, ie n ; q)G(k, ie n )G(k + q, ie n ). (B-6) 

k,?e n 

Appendix C: Quartic AF mode-mode coupling strength 

In this appendix, the quartic coupling strength among the four AF spin fluctuation 
modes is calculated. This coupling strength is equivalent to the fourth-order coefficient of 
the Ginzburg-Landau free energy for an AF order parameter. 

We assume the same single-particle dispersion used in Appendix A. The quartic coupling 
74 of AF spin fluctuation modes at Q* = (tt, it) is evaluated as 

1 



74 



k n (ie„ - e(k) + ^i) 2 {ie n - e(k + Q*) + (j,) 2 

TY I 

^ (ie„ - e(k) + ^) 2 {ie n + e(k) + n) 2 

roc i 
T \ Ap / de 



/c 



(ie n - e + n) 2 (ie n + e + fx) 



2 



n>0 



1 1 

+ 



(ie n +fi) 3 (ie n -n) 3 _ 



^ r R e ^) fl + M), (C-l) 



(4vrT) 2 r V 2 2 7rT, 
where V^ 2 ^ is the tetra-gamma function. Re^( 2 )(| + ix) changes its sign at x ~ 0.304; i.e., 
Reip^d + ix) < (> 0) for x < 0.304 (x > 0.304). 56 Thus, 74 changes its sign at \n\ ~ 1.91T 
from positive to negative as increases; i.e., 74 is expected to be negative in the low- 
temperature region T < 0.524|^| in the underdoped region of the simple Hubbard model with 
t' = 0. The negativity of the quartic coupling constant 74 leads to a first-order phase transition 
and makes higher-order AL terms positive. 
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Note added in proof: 

After the present paper had been accepted for publication, we realized that we had missed 
to refer previous works which reported an experimental aspects related to the physics discussed 
in the present paper. 

Effects of defects in correlated superconductors are reviewed in H.Alloul, J. Bobroff, M. 
Gabay, and P. Hirschfeld: Rev. Mod. Phys. 81 (2009) 45, in which the role of disorder in 
single-layered curate Hgl201 compound is discussed. The importance of taking into account 
the disorder for understanding the phase diagram, especially apparent competition between 
magnetism and superconductivity, of double-layered cuprate YBCO systems was discussed 
from experimental side in F. Rullier-Albenque, H. Alloul, F. Balakirev, and C. Proust: Euro- 
phys. Lett. 81 (2008) 37008. 

We thank Prof. Henri Alloul for reminding us of those works. 
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